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The parametric formulas are obtained whereby the refractive index and the light group 
velocity could be calculated for any water state (salinity, pressure, temperature) inside of 
the light wavelength region from 200nm to 1200nm. These formulas are suitable to be used 
in underwater neutrino telescope activity. 

Introduction 

In underwater Cherenkov neutrino telescopes /l-j-5/ water acts as a radiator of 
Cherenkov photons generated by charge particles. Cherenkov photons are emitted 
in cone with half-angle defined by the mass and momentum of the particle and the 
refractive index of the medium (water in our case) 

cosQ = l/(/3 • n(X)), (1) 

where (3 is velocity of particle, n(X) is refractive index and A is wavelength. 

Because the medium (water) is disperseve, the group velocity (velocity of light 
signal transmission in transparent material) is different from phase velocity 
Vph(X) = c/n(X) and being characterised by a group index of this medium 

v gr (X) = c Av(A), (2) 

n gr (X) = n(X) — X ■ (dn/dX), (3) 

where c is the light velocity in vacuum. 

This reality is taken into account for example in construction of a large-area Time 
Of Flight(TOF) counters. The group index of scintillator is n gr = 1.75 compared 
with refractive index n = 1.5 /6/. In the DUMAND-type installation activity a 
phase velocity was used for estimation of a photon TOF till it was demonstrated 
jl I that such approach leads to errors in results of TOF-calculations. For photon 
passed through 100m water layer these errors vary from 10ns to 18ns when the light 
wavelength varies from 500nm to 350nm. A disregard of the v gr — v p h difference does 
not affect errors in track reconstruction for the existing projects /8/, but, as Baikal 
Neutrino Telescope experience shows, only v gr must be used for time-calibration of 
the detector modules with the help of an outside laser and its value have to be known 
with a quite high accuracy. 

The measurement of the group velocity may be performed with a help of detector 
modules lighted by an outside laser which shines through water /9/. Comparison 
of the experimental result with the value of this physical constant gives an unique 
possibility to verify the accuracy and reliability of all time measuring units and all 



time calibration procedures. To prove the height level of the TOF measuring the 
accuracy of v gr calculation have to be higher than this level and therefore all factors 
which affect the n gr value (wavelength dependence, temperature, pressure, salinity) 
must be taken into account. 

1. Dispersion formulas for pure and sea water 

In applied optics a calculation of a refraction index for any wavelength in the 
restricted wavelength region is performed using so called dispersion formulas. These 
formulas have the same structure for different class of material but differ usually in the 
formula coefficients. Thus for example for optical colourless glass in the wavelength 
interval 365H-10139 nm the following formula is used /10/: 

n 2 = A, + A 2 ■ A 2 + A 3 • A" 2 + A 4 • A" 4 + A 5 • A~ 6 + A 6 ■ X' 8 , (4) 

where are the coefficients of the dispersion formula, which depend on the glass 
class and each class is characterised by one's own set of Ai /10/. 

There is abundance of table data/10-12/ which contain the results of the refrac- 
tion index measurements performed for pure (salinity S=0 °/oo) water, atmospheric 
pressure P = P a tm and temperature t = 20°C. The formula (5) represent the best fit 
to these data and the lower solid line 1 on Fig.l shows the result of the refraction 
index calculation by means of this expression. The points on this line correspond to 
table data used in the fitting procedure. 

n 2 (X) = 1.7527 - 2.55 • 1(T 3 • A 4 - 5.0 • 10~ 3 • A 2 + 9.9 • 10~ 3 • A~ 2 - 

-4.0 • 10" 4 • A" 4 + 2.9 • 10~ 5 • A~ 6 - 5.0 • 10~ 7 • A" 8 , (5) 

where A in mkm and index '0' means distilled water with t=20°C and P=Patm. By 
comparison with expression (4) the term A ■ A 4 was added to extend the region in 
which dispersion formula is correct since A=214nm up to A=1256nm. 

Salinity, pressure and temperature affect the value of the water refraction index. 
The sets of the same black points which are also shown in Fig.l correspond to the 
sets of the table data /12/ for different states (S, P, t) of sea water. The dependencies 
n versus A for these sets are described by functions 

n(X, S, P, t) = n (X, S , P , t ) + An(S, P, t) (6) 

with quite high precision. The last fact permits to conclude, that in distinct with 
an optical glass, the dispersion formula coefficients are the same for distilled water, 
sea water and compressed water and they do not depend on temperature. Thus to 
calculate a refraction index of any water for any wavelength the dispersion formula 
(5) may be used as an universal one and just S, P, t -corrections (An(S, P, t j) must 
be taken into account which do not depend on wavelength. 
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Table data for pure water (t = 20°C,P = P atm ). 
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Table data for sea water(S = 35 °/ 00 , t = 5°C) /1 2/ 
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Phc. 1: Dispersion formulas and table data. The lines (1-5) were calculated by means of the 
dispersion formula(5) and line 6 by means polinomial /14/. 1 : n=no; 2 : n=no+0.0081; 
3 : n=n +0.0132; 4 : n=n +0.0182; 5 : n=n +0.0266. 



2. S, P, t -corrections of water refractive index 

In the simplest approach /7, 8/ the Lorentz-Lorentz formula may be used to get such 
corrections: 

n 2 ~l 

n 2 + 2'p 3'M' ae ' {) 
where M is the molecular weight, Na - Avogadro number, a e - polarizability and 

n 2 - 1 1 



R 



P ~ 2 



n 2 + 2 p 



(8) 
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is the specific refractive index. If we suppose that the right part of the expression (7) 
is constant, the change in refractive index with density is given by 

An (n 2 + 2)-(n 2 -l) 

Ap 6-n-p {) 

and calculation of Ap(S, P, t) is enough for the estimation of An(S, P, t). 

Analysis of table data/ 12/ shows some dependence R p versus t and versus S 
and so the described approach becomes not quite correct. Direct comparison of the 
measured increasing A n exp which occurs due to water temperature fall from 20°C 
to t°C with the value A n^_£ (formula (7)) shows that the measured results exceed 
the estimated values on 18%, 23%, 29% and 37% for temperatures 15°C, 10°C, 5°C 
and 0°C correspondingly. Therefore in the given calculations the experimental results 
are used as fully as possible and the formula (9) may be used, if such results are not 
available, but only for An p correction, when density is changed directly due to the 
pressure increasing. To calculate the corrections it is supposed that 

An(S, P, t) = An St (S, P , t) + An p (S, P, t). (10) 



The refractive index table data which have been measured for sea water /12/ were 
used to get some mathematical expressions to describe the dependencies n versus S, t 
and P. The points on Fig. 2 correspond to sea water refractive index values measured 
under atmospheric pressure for photons with A = 583. 9nm /12/ and the lines are 
calculated using the expression (11), which is the best fit to these points: 

n(S,P ,t) = 1.3340 + 0.19705- 1(T 3 • S- 

-(0.15655 + 0.010045 • S - 0.41077 • 10~ 4 • S 2 ) ■ 10~ 4 • t- 

-(0.17585 - 0.14300 • 10" 2 • S + 0.12525 • 10" 4 • S 2 ) ■ 10~ 5 • t 2 . (11) 

Because of the salinity varies inside the region (35 ±4)°/ 00 in projects planned to 
be installed in a deep sea, the table data /12/ are used which contain refractive index 
values (A = 501. Inm) measured under high pressure in water with S = 35°/oo- These 
points for which the sea depth values (h(m)) are used instead the pressure values 
are shown on Fig. 3 and dependence of the refractive index versus h and versus t is 
described well by the expression (12) : 

n(35, h, t) = 1.3448 + 0.15024 • HT 5 • h - 0.13122 • 10~ 10 • h 2 - 

-(0.50560 + 0.80127 • 10' 4 • h - 0.99739 • 10~ 9 • h 2 ) ■ HT 4 • t- 
-(1.3118 - 0.11469 • 10~ 3 • h + 0.18712 • 10~ 8 • h 2 ) ■ 10~ 6 • t 2 . (12) 

The expressions (11) and (12) may be used to calculate (S, t)-correction and P^- 
correction correspondingly : 

An St = n(S,Po,t)-n{0,Po,to), (13) 
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An p = n(35, Ph, t) — n(35, Pq, t). 



(14) 



A pure water compressibility is higher than the one for salt water and, therefore, 
according (9) one has An p (0,P, t) > An p (S, P,t). To account it (in the absence of 
direct experimental data) the depth correction for pure water may be calculated as 

An p (0, P h , t) = An p (35, P h , t) ■ K s , (15) 




Phc. 2: Refractive index (A = 589. 3nm) of sea water under atmospheric pressure /12/. 



■— ■ - ~-> yj — • 
















1.3625 






o 


t=0.03°C 








1.36 






□ 


t= 1 5.02 ° C 






. - A' " 


1 .3575 








t=29.98 ° C 








1.355 












^ ' ' ' ' - A " 




1.3525 
















1.35 
















1 .3475 








. , A' ' ' ' 








1.345 












S 3 5 / oo 




1 .3425 
















1.34 






I , 










I 







2000 


4000 6000 




8000 1 0OOO 1 2000 


14000 



depth m 



Phc. 3: Refractive index (A = 501. 7nm) of sea water (S=35°/oo) for various depths /12/. 

where K s = Ap(0, P h ,t)/ A p(35,P^,t). Analysis of the table data /12/ gives 
X s = 1-07 ± 0.01 for t = (0 + 5)°C and /i = (0 14000)m. 
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3. An(S, P, t) -corrections and group light velocity for existing projects 

The refractive index corrections caused by environmental parameters (S,Ph,t) of 
detector sites, and values of the light group velocity for A=475 nm which have been 
calculated by means expressions 5-j-6 and 10-i-14 are put in the table. The points on 
Fig. 4 reflect the results of such calculations for the various values of light wavelength. 



Detector 


s (%») 


h (m) 


t(°c) 


An(S,P h ,t) 


v gr (10 10 cm/s) 


DUMAND 


35 


4600 


1.3 


0.0147 


2.1729 


NESTOR 


38.7 


3800 


14 


0.0129 


2.1757 


ANTARES 


33.65 


2500 


13.2 


0.0103 


2.1799 


BAIKAL 


0.12 


1150 


3.4 


0.0027 


2.1919 




Phc. 4: The dependencies of calculated light velocity versus light wavelength for the existing 
underwater neutrino telescope projects. The points and upper line have been calculated by 
means the expressions (2-=-3), (5-r6) and (10-r-14). The stars correspond to the experimental 
results/9,14/. The other lines are calculated by means the formula(16). 

It seems to be completely justified to simplify the calculation procedure by using 
the following formula which is the best fit for these points: 

v gr (X, S, P, t) = c/(1.333 + 0.57 • 10~ 2 • A~ 2 + 0.3 • 10" 3 • A~ 4 + An(S, P, t)). (16) 
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4. Comparison with the experimental result 



The first direct measurements of the light (A=475nm) signal velocity (V) have been 
performed in the water volume of Baikal Neutrino Telescope (NT-200) /9/. These 
measurements did not require any add devices and any changes in NT-200 standard 
regime. Data of usual laser calibration run has been used to do this. During such 
runs NT-200 Channels (Ch) detect the flashes of laser, which has attached to NT- 
200 central string on 20m lower than the lowest Ch. Laser generates five flash series 
with the decreased photon intensity. Only central string data have been taken into 
account to avoid the influence of the neighbour strings co-ordinates ambiguities. Mean 
measured values of Time Difference (T y -) for any of two Ch's (45 pairs in the total) and 
the (measured before) spaces between Ch's attachment points and laser attachment 
point were used as an input information. Increasing of Ch time delay with the laser 
flash intensity decreasing (time walk effect) was observed for all Ch's. The special 
procedure /9/ was applied to compensate an influence of this effect on the final result. 
Analysis of the ^-distributions measured for different SV, = Si — Sj (Si is the distance 
from laser to Chj) has demonstrated no dependence of mean velocity values < V > 
on Sij (Fig 5A). For selection criterion Si — Sj > 12.5 m the expected value of the 
laser signal velocity 

V = (2.178 ± 0.014) • 10 10 cm/s 

had been obtained using data of 30 pairs of Ch's which satisfy this criterion (Fig 5B). 
This result is comparable with predicted value 2.1919-10 10 cm/s (see table and Fig. 4), 
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Phc. 5: The results of in situ measurements of light signal velocity /9/. 
the distinction does not exceed one standard deviation (~ 0.7% of the measured value) 
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and it is possible to conclude with high confidence level that the experimental light 
signal velocity value agrees with the calculated value of light group velocity and does 
not agree with the phase velocity value. To verify an influence of salinity, pressure 
and temperature on velocity the precision of measurements has to be increased. 

Next more extensive and more precise measurements have been done with the 
experimental string developed near NT-200 site to test new deep underwater tech- 
niques /14/. Light Beacons supplied short flashes with fixed photon intensity were 
used as light sources in these measurements. UV(A =370 ± 6 nm),blue (A =470 ± 
11 nm) and green (A =520 ± 17 nm) LEDs were used in Light Beacon. LED flashes 
were detected by four standard NT-200 Channels attached to experimental string. 
The measurements of the time intervals between trigger signals of LED drivers and 
time responses of Chennels, the time delay in electronics and cables (obtained in the 
laboratory) and the distances between light source and Channels (31m, 34m, 58m, 
64m) were used to get the light signal velocity. In contrast with first experiment 
there were no ability to estimate the influence of the time walk effect, nevertheless 
the results : (2.148±0.010) • 10 10 cm/s (A = 370 ± 6 nm); (2.193 ± 0.009) • 10 10 cm/s 
(A = 470±11 nm); (2.206±0.009) • 10 10 cm/s (A = 520± 0.017 nm) are in the excellent 
agreement with the predicted values which are 2.1516-10 10 cm/s, 2.1920-10 10 cm/s 
and 2.2029-10 10 cm/s correspondingly (see also Fig4.). 

5. Summary and conclusion 

The measured properties /10-j-12/ of the water refractive index as functions of 
light wavelength(A), salinity(S), pressure(P) and temperature(t) have been employed 
in the performed analysis as fully as possible. 

The structure of dispersion formula for water was applied which is similar with 
dispersion formula structure for an optical glass /10 /. This structure was used to 
describe sucsesfully (see formula (5) and Fig.l) the experimental refractive index 
dependence on the light wavelength in the interval (200-=-1250)nm for distilled water 
with t=20°C and P=P aim /10-=-12/. The polynomial using/14/ gives true result only 
in the interval (300-^600)nm (see Fig.l) in which Cherenkov light contributes to the 
signal of a muon. 

Analysis of the extensive table data jYlj shows that formula (5) may be also 
successfully used (see fig. 1) for water with other properties (S, P, t) just by means 
adding some correction Anspt{S, P, t) which does not depend on A. 

The assumption is used that equality Ansp^An^+Anp is sufficiently true in 
the considered case, and these corrections are computed for any water states using 
expressions (ll-=-15). 

To simplify calculations of the light group velocity the suitable formula (16) has 
been selected which quite satisfactorily describe the computed group velocity values 
in the light wavelength interval 200-^1250 nm (see Fig. 4). 

The comparison of the group velocity value calculated for Baikal Neutrino Tele- 
scope site with the results of in situ measurement /9, 14/ demonstrates good agree- 
ment, the distinctions do not exceed one standard deviation (~ 0.7% of measured 
value) /9/ and 0.4 standard deviation /14/ while the difference between calculat- 
ed phase velocity value and the measured light velocity amounts to four standard 
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deviations. Moreover the results /14/ permit to confirm a validity of Av 9r (S, P, t) 
estimations (see Fig. 4). 

The obtained parametric expressions (5), (6), (10-=-16) could be used to calculate 
a velocity of monochromatic light signal for calibration of TOF system by means of 
a power laser lighting through water. These expressions are also suitable for using in 
simulation of events detection and in reconstruction procedure for existed Underwater 
Neutrino Telescope projects and for future lkm 3 Neutrino Telescopes /13/. 
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